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FREE-VIBRATION CHARACTERISTICS OF A LARGE SPLIT-BLANKET 
SOLAR ARRAY IN A 1-g FIELD 
by Francis J. Shaker 
Lewis Research Center 

SUMMARY 

Two methods for studying the free -vibration characteristics of a large split -blanket 
solar array in both a 0-g and a 1-g cantilevered configuration are presented. The 0-g 
configuration corresponds to an in-orbit configuration of the array; the 1-g configuration 
is a typical ground-test configuration. The first method applies the equations of con- 
tinuum mechanics to determine the mode shapes and frequencies of the array; the sec- 
ond method uses the Rayleigh -Ritz approach. In the Rayleigh -Ritz method the array 
displacements are represented by string modes and cantilevered beam modes. The re- 
sults of this investigation are summarized by a series of graphs illustrating the effects 
of various array parameters on the mode shapes and frequencies of the system. The 
results of the two methods are also compared in tabvilar form. 


INTRODUCTION 

To predict the dynamic behavior of a spacecraft with large, flexible solar arrays, 
it is first necessary to describe the structviral dynamic characteristics of the array . 
These characteristics can be determined once the cantilevered modes and frequencies 
of the flexible solar array in a 0-g environment are known (refs. 1 and 2) . 

In addition to the 0-g modes and frequencies, it sometimes becomes necessary to 
predict the dynamic behavior of the array while it is suspended vertically in a 1-g field. 
This is the usual ground-test configuration for studying the dynamic behavior of large 
solar arrays (refs . 3 and 4) . In this configuration the blanket weight has a significant 
effect on the dynamics of the array and must be included when predicting the behavior of 
the system. 

The cantilevered modes and frequencies of a split -blanket solar array have been 
studied by several investigators in recent years. As shown in figiire 1 a split -blanket 
solar array consists of a central support boom, a two-piece blanket substrate with solar 



cells mounted on one side, and a leading-edge member that transfers a tension load to 
the blanket from the boom. For this lype of array the boom centerline lies in the plane 
of the blanket. In reference 5 the cantilevered modes and frequencies of a split -blanket 
array in a 0-g field were investigated by solving the differential equations governing the 
motion. This method results in transcendental equations that can be solved numerically 
for the frequencies. In reference 6 the finite element method was used to find the canti- 
levered modes and frequencies of this type of array. Beam -column elements were used 
to describe the array boom, and membrane elements to describe the blanket. Seventeen 
and 26 degrees of freedom were used to determine the fundamental antisymmetric (tor- 
sional) and symmetric (bending) modes and frequencies, respectively. Using the finite 
element method to determine higher modes and frequencies would require a large num- 
ber of degrees of freedom to obtain the same accuracy. In reference 7 the modes and 
frequencies of a solar array in which the blanket is offset from the centerline of the 
boom are investigated by using a Rayleigh -Ritz approach. The assumed modes used to 
construct the Rayleigh -Ritz solutions were polynomials to the fourth order. 

The purpose of the present report is twofold. First, it presents the exact solution 
for the modes and frequencies of a split-blanket solar array in a 1-g field and compares 
this solution with a Rayleigh-Ritz approximate solution. Second, by using these solu- 
tions it investigates the effects of various solar array parameters (mass, blanket ten- 
sion, boom stiffness, etc.) on the modes and frequencies of large solar arrays in both 
a 0-g and a 1-g configuration. In developing the Rayleigh-Ritz solution, cantilevered 
beam modes and string modes are used to represent the boom and blanket displacements, 
respectively. Using these modes will ensure convergence for a relatively small number 
of modes or degrees of freedom. The effects of solar array parameters on the modes 
and frequencies are illustrated by a series of graphs, and comparisons between exact 
and approximate solutions are shown in tabular form. 


THEORETICAL DEVELOPMENT OF NORMAL MODES AND FREQUENCIES 
OF LARGE SPLIT -BLANKET SOLAR ARRAY 

Continuum Mechanics Approach 

Basic assumptions . - For purposes of analysis a large split -blanket solar array is 
idealized, as shown in figure 1. This figure shows the array consisting of three compo- 
nents; a center boom that supports the array; a membrane substrate with solar cells 
attached to one side (hereinafter referred to as the blanket) ; and a bar at the tip of the 
boom that transfers a tension load P from the boom to the substrate. The displace- 
ments of the boom and blanket, normal to the plane of the blanket, are denoted by V(x,t) 
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and W(x,y,t), respectively. In developing the equations of motion for the array the fol- 
lowing assumptions were made: 

(1) The bending stiffness of the blanket, normal to its plane, is negligible so that 
the blanket behaves like a membrane in this direction. 

(2) The tension distribution is uniform across the width of the blanket (i.e. , the 
tip piece is perfectly rigid) . 

(3) Displacements are small, so that small -displacement theory is valid. 

(4) Boom weight is negligible and the shear center coincides with the neutral axis of 
the boom. 

Based on these assumptions the equations describing the motion of the array were de- 
veloped. (All symbols are defined in appendix A.) 

Equations of motion . - The forces acting on an element of the blanket displaced an 
amount W(x,y,t) from its static equilibriuip configuration are shown in figure 2(a) . 
Applying Newton's second law of motion to this element yields the following equation: 

T 9 /9W\ 

3x \ 3x / 9x dx b 

where b is the blanket width and is the mass per unit length of blanket . It can 
readily be seen that an alternative form of equation (1) is given by 

-L (t ^ (2) 

8x V 0x/ b g^2 

Now for a blanket hanging vertically the tension at any point x will be a superposition of 
the preload P transferred to the blanket from the boom and the weight of the blanket 
below the point x. That is. 


T(x) 


P / ^ _ ^m^ 

b \ P P 


(3) 


where is the total blanket weight. In view of equation (3), equation (1) is trans- 

formed by making the following change of variables : 




( 4 ) 


From equations (1) and (4) then 
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(5) 


a^2 ? a? ■ P 3^2 


Equation (5) represents the desired form for the equation of motion of the blanket. 

The equation of motion of the boom can be developed by applying Newton's second 
law of motion to the beam element shown in figure 2(b) and by using the following force - 
displacement relations from elementary beam theory: 



El ax 


(6a) 


M (x) = El 

y 2 

ax*^ 


(6b) 


In this manner it can be shown (ref. 8) that the equation governing the bending motion of 
the boom is given by 



P a^v 

El ax2 


Pb a^v_ 

El 3^2 " 


(7) 


In addition to the bending motion described by equation (7) , the boom can also experience 
a rotational motion about its centerline . The equation governing this motion is developed 
in numerous texts on vibration theory (e.g. , ref. 9) and is given by 

(8) 

ax^ at^ 

where 9 is the rotational angle of the boom cross section, I^ is the mass polar moment 
of inertia per unit length, and JG is the torsional stiffness of the boom. The relation 
between the torsional moment M and the angular rotation 9 used in developing equa- 
tion (8) is 


M^(x) = JG ^ (9) 

^ ax 

Equations (6) to (9) represent the required relations for the boom. 
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The final set of equations are the equations of motion for the tip piece. The forces 
acting on the tip piece are shown in figure 3. Applying Newton's second law of motion 
for forces in the z-direction and moments about an axis parallel to the x-axis and pass- 
ing through the center of gravity of the tip piece yields the following two equations; 





p 3W(i,y,t) ^ , . 3^V(i,t) 

» 3* 3t2 


= 0 


M^(^t) + 



awq.y.t) 

9x 


dx + I, = 0 

IP n 

dt^ 


( 10 ) 


( 11 ) 


where is the mass of the tip piece and I|.p is the mass moment of inertia about its 

center of gravity. Eliminating the force and moment in equations (10) and (11) by using 
equations (6a) and (9) and changing the variables in W from x to ? by using equa- 
tion (4) transforms equations (10) and (11) into the following equations: 


8^v(3,t) , p av(^.t) , 

3x3 El ax 



M-i.y^dy 

a? 


El at2 


( 12 ) 


ae(i ,t) . 1 /P + (1 

ax 2^b\JG/\ p / / 9? JG ^^.2 

t/-b/2 

Equations (12) and (13) represent the final form of the equations of motion of the tip 
piece. These equations together with equations (5), (7), and (8) represent the motion 
equations for the complete solar array. The displacement variables in these equations 
must satisfy certain boundary and compatibility relations. These conditions are given 
next. 

Boundary conditions and compatibility relations. - At the fixed end of the array the 
displacements and rotations of the array elements are all zero. At this end, x = 0 and 
it follows from equation (4) that ^ = yi + (W^^^/P) . Thus, the boundary conditions 

at the fixed end will be as follows : 
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V(0,t) =0 


>1 


— (0,t) - 0 
dX 


> 


(14) 


0(0, t) -0 


w(?Q,y,t) 


At the free end of the array the displacements and rotations of the components must be 
compatible. At this end x = I ; and from equation (4), ? = 1. In addition, the moment 
at the tip of the boom (given by eq. (6b)) is zero. Thus, the boundary and compatibility 
relations at the free end are as follows: 


0{l,t) a 
W(l,0,t) = V(^,t) 
W(l,y,t) = W(1,0, t) + ya 

dx^ -> 


(15) 


Equations (12) to (15) represent the complete set of relations that must be satisfied by 
the solutions to equations (5), (7), and (8). 

Separation of variables . - The first step in the solution of the motion equations is 
to eliminate the y and t variables by assuming a harmonic solution of the following 
form: 


V(x, t) = v(x)e^‘^^ 
W(?,y,t) = [w(?) +y<p(0y‘^* > 
0(x,t) = 0(x)e^“^ ^ 


(16) 


From equations (16) and (5), (7), (8), and (12) to (15), two independent sets of equations 
and corresponding boundary conditions are obtained. These equations can be nondimen- 
sionalized by making the following substitutions: 
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X = ^, v = l, w = ^, P = ^, = 

I I I * P 

m 


El ^ 


, n-i. - 

El ^ JG 




_ M.„ - H - - Ih 

M. = -i£, b = ^, I. = — , I. = 

Mj, ^ i/r V.2 ^ k2 

^ Pm^ 


2 

z m 


«! = 


k2 


(17) 


where is the total mass of the blanket and Mj^ is the total mass of the beam. The 

first nondimensional set of equations is designated as the bending equations and is given 
by 


! « 


dfv.j4-,0 
dx^ dx^ 


3!vm . -^2 dv^ ^ 1 ^ dwjli ^ 


dx^ 


dx 2 P d? 


with boundary conditions 


(18a) 


(18b) 


(19) 


o 

11 

o' 

l"> 

(20a) 

dv(0) _ 0 

(20b) 

dx 


w(?q) = 0 

(20c) 
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w(l) = v(l) 


d^v(l) - 0 

dx2 


The second set of equations is called the torsion equations and is given by 




d? 


2 ? d? 


^.Ibk2^0 = O 


dx‘ 


k? 


. J_4. Mli . r kf^0(i) = 0 

dx 24 P d? tp t t 


(21) 


(22a) 


(22b) 


(23) 


with boundary conditions 


0(0) = 0 

(24a) 

il 

O 

(24b) 

0(1) = <p{i) 

(24c) 


Equations (18) and (19) and boundary conditions (20) and (21) define the bending or sym- 
metric motions of the array in which the blanket displacements are independent of y. 
This type of displacement field is illustrated in figure 4(a). Similarly, equations (22) 
and (23) and boundary conditions (24) define the torsional or antisymmetric motions of 
the array, in which the blanket displacements vary linearly with y. Figure 4(b) illus- 
trates this type of displacement field. 

Solution to bending equat ions. - To obtain a solution to the bending equations, a 
solution to equations (18) that satisfies equations (19) to (21) must be determined. Now, 
equation (18a) is Bessel's equation of zero order, and its solution satisfying boundary 
condition (20c) is (ref. 10) 

w(0 = c[Yq(2/Cj^2^0)Jq(2kiF^O - Jo(2«j^2^q)Yq(2/<^^2^)] (25) 
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where Jq and Yq are Bessel functions of the first and second kind, respectively, of 
zero order and C is an unknown constant. Also, it can be shown by direct substitution 
that the solution to equation (18b), satisfying boundary conditions (20a) and (20b), is 

v(x) = A(q ’2 cosh ojjX - o/g cos Qfgx) + B(ag sinh a^x - a-^ sin (26) 

where 



and A and B are unknown constants. Now equations (25) and (26) contain three un- 
known constants that can be evaluated from the three remaining conditions given by 
equations (19) and (21). Utilizing equations (25) and (26) changes equations (21) to 

A(q !2 cosh Q!j - o >2 cos cug) + ^ioi2 “l ~ “l “2^ “ CFj(j3) = 0 (28) 

A^cvj cosh Qtj + cvg cos sinh otj + * 2 ) ~ ^ 


where 




(30) 


The third equation is obtained from equation (19) by using equations (25) and (26). After 
some manipulation with Bessel functions, this eqxiation can be written as 

A sinh Qfj - o; 2 a '2 ^tp^^“2 “l " “2 cos 02 )] 

+ B cosh o!j^ + cos + ^tp^^“2 ^1 ” ‘^l *^2il 


fi? -1 

: ^ fj/3) = 0 


(31) 
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where FgC/S) is defined in terms of Bessel functions of the zero and first orders as 



( 32 ) 


Equations (28), (29), and (31) represent three equations in three unknowns. A, B, and 
and C, For a nontrivial solution the determinant of the coefficients of the unknowns 
must be zero. Expanding this determinant and simplifying yield the following equation: 


FjWMj/.SfjW [{4*aiy 


(q!o sinh a, cos 


01-. 


cosh sin ofg 


>] 


+ F 


l(/3)[2/S^ 


+ 


/3^(2/3^ 


+ k^)cosh 


a-i cos ac 




sinh sin a 




(33) 


Equation (33) is the characteristic equation for the bending vibrational frequencies of a 
split-blanket solar array in a 1-g field. When the mass parameters and Mj^ and 

the load parameters P and k are specified, the frequency parameters |3 can be de- 
termined numerically. Note that equation (33) is valid, provided the boom does not 
buckle. That is, the load P must be less than the critical buckling load of the array. 
Development of this critical buckling load is presented in appendix A. For P equal to 
Per, it can be shown that equation (33) yields = 0. 

Solution to torsion equations . - The torsion equations can be solved in a manner 
similar to that used for the bending equations. It can be verified by direct substitution 
that the solution to equations (22) satisfying the boundary conditions (24a) and (24b) is 
given by 


0 = Aj sin(l^/2j^^/3tx) 

<p = Bi[Yo(2/C2?o)Jo( 2«2?) - Jo(2/C2?o)Yo(2«2?)] 


(34) 


where «2 = The constants Aj and B^ are such that conditions given by equa- 

tions (23) and (24c) must be satisfied. By using equations (34) these equations become 


Aj sin^jj^^j./3^.) - BjF2(/3j.) - 0 


k? 


(35) 
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where 


Fjdt) = Y„(2«2?„)J„(2«2) - J(,(2r2C„)Y„(2r2) 


” J (L ) 
0 2 




(36) 


Equating the determinant of the coefficients of the unknowns Aj and Bj in equa- 
tions (35) to zero yields 


12k, 


F ^0^)sin{ly^ 




T1/2 
— cos 


(ifa^Vt) - Vt Vt) 


= 0 


(37) 


Equation (37) is the characteristic equation for the torsional frequencies of the array in 
a 1-g field. Once the inertia parameters and 1^^ and the load parameters k^. and 
P are specified, the frequency parameters j3^ can be found by solving 
f the mass of the boom is negligible such that sin^I^^^ 

» 1 , then equation (37) simplifies somewhat to 


numerically . I: 
and cos^I^'^^j, 


this equation 


'^2 



= 0 


(38) 


Equation (38) is applicable to current large solar array designs. 

Degenerat e case - 0-g configimation . - For the case of a solar array in a 0-g field 
(i.e. , an in-orbit configuration) the characteristic equations for the bending and tor- 
sional frequencies can be determined from equations (33) and (38) by taking the limit of 
the functions in these equations as the blanket weight approaches zero. It is shown 

in appendix B that as approaches 0 the following relations hold: 


Fl(/3) = 




sin o! . 


(39a) 


Fa(^ = 


z2).l/2 


COS Of, 


(39b) 
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1 


(40a) 






^4(%) = cos 

^ L 1/0 


7T? 


1/2 


(40b) 


where 




“3 = 


2 -^ 
k 


(41) 


Substituting equations (39) into equation (33) yields the following characteristic equation 
for the bending frequencies of the array in a 0-g configuration: 

(^tp^^ sin Qfg - k^ag cos + ^^^(ag sinh cxj cos a 2 - aj cosh oij sin 012 )^ 


+ sin dig [ 2 ^ + /3^(2^ + k^cosh cos a 2 - sinh sin 02 ] = 0 


(42) 


Similarly, the characteristic equation for torsion is obtained by substituting equa- 
tions (40) into equation (38), which yields 


cos /3t + 12 / — - a sin /3t = 0 

m 


(43) 


It can be shown that equations (42) and (43) are analogous to the frequency equations in 
reference 3, which were developed in a straightforward manner for the 0-g configura- 
tion. 

Mode shapes and orthogonality conditi ons. - To determine the bending mode shapes 
of the solar array, assume that the frequency parameter for the n^^ mode is known. 
Corresponding to there will be an and oi2n‘ Starting with equations (25) 

and (26) the boom and blanket displacements in the mode will be 


Vn(x) = A^(« 2 n cosh x - cos o< 2 nX) + B^(Q- 2 n sinh a^x - sin o, 2 nX) 


(44) 
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(45) 




The mode shapes can be determined from equations (44) and (45) by expressing and 
Cjj in terms of using the boundary conditions expressed by equations (28) and (29) . 
From these equations then 


B ~(^ ^n °‘ln ^ ^2n °^2n) ^ 

“fn “in + “ln“2n “2n 


Cn = 




Q! 2 n(cosh - cos 


(46) 


(“L “In + “L “2n)“2 n “in ' “in “2n^ 

+ “ln“2n “2n 


(47) 


“In “In 


If equations (46) and (47) are substituted into equations (44) and (45), the n^*^ mode shape 
for an array in a 1-g field can be determined to within an arbitrary factor A^. For the 
0-g field the boom equation given by equations (44) and (46) remains unchanged. The 
blanket equation for this case becomes 


^n^^^ " ^n “3n^ 

In addition, for the 0-g field the boundary condition given by equation (28) yields the fol- 
lowing result for 


D 


n 


sin 0 - 3 ^ 


(49) 


J.L. 

Thus, equations (44), (46), (48), and (49) determine the mode for an array in a 0-g 
field. The mode shapes for the torsional cases can be developed in an analogous man- 
ner. 

The orthogonality relations for the solar array can be developed in a straightforward 
manner. Because of the length of this development it is given in appendix C. 

Summary of continuum mec hanics approach . - In the previous sections the mode 
shapes and frequencies of a split -blanket solar array were determined by solving the 
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differential equations governing the motion of the system. The results showed that the 
array will exhibit symmetric or bending modes and antisymmetric or torsional modes of 
vibration. The frequencies of the bending modes can be found by solving equation (33) 
or (42) for an array in a 1-g or a 0-g field, respectively. Similarly, the torsional mode 
frequencies can be determined from equation (37), (38), or (43) depending on the partic- 
ular case of interest. In all of these cases a highly transcendental equation involving 
Bessel functions and/or trigonometric and hyperbolic functions must be solved numer- 
ically to obtain a solution. In the following section an alternative approach is developed 
that uses the Rayleigh -Ritz method. This approach can readily be extended to include 
more complex arrays for which exact solutions cannot be found . 


Rayleigh -Ritz Approach 

In the Rayleigh -Ritz method (ref. 9) the strain energy and kinetic energy of the sys- 
tem are expressed in terms of (1) assumed modes that are functions of the spatial co- 
ordinates and (2) unknown generalized coordinates that are functions of time. By using 
Lagrange's equation and these energy expressions a set of second-order, linear, differ- 
ential equations in terms of the generalized coordinates can be obtained. This set of 
equations is then reduced to an eigenvalue problem from which the modes and frequencies 
can be determined. 

The success of the method is highly dependent on choosing a proper finite set of as- 
sumed modes that will ensure accurate results. The method, as applied to a large solar 
array, is now developed for the array bending vibrations and the results compared with 
the exact results of the previous section. 

Strain energy of system in bending. - Let be the uniform stress resultant in 

the blanket when it is in its static equilibrium configuration. During the motion of the 
blanket about this equilibrium configuration the resultant stress will change to As- 
sume that the stress induced by motion is small relative to so that N„ . 

Under these conditions it can be shown (ref. 11) that for small displacements the change 
in strain energy in the blanket is given by 

ph/2 pi 2 

J-h/2 Jq 

For an array hanging vertically in a 1-g field the resultant stress is given by equa- 
tion (3) as 
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n(o) = £ a + 1 . 

* b V P P/ 


From equations (50) and (51) the blanket strain energy in nondimensional form is 


= 

m 


/ V P P/Vax/ 


where W = W/l. Similarly, it can be shown (ref. 12) that the strain energy in the boom 
is given by 


f 

^21 j 


cb?-k2 


/ 


where V =W /I . The total strain energy in the array is then 


U = + U. 

m D 


Kine tic energy of system in bending . - The kinetic energy of the array consists of 
the kinetic energy of the boom, blanket, and tip piece. The kinetic energy of the blanket 


T = 
m 


1 ^ fdW 

2 / ^“(at. 


dx = 

2 “ 


/ 


aw\^ 


Similarly, the kinetic energy of the boom Tj^ is 




dx 


Lastly, the kinetic energy of the tip piece T^p is 
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(57) 



tp" 



2 


The total kinetic energy of the array is then given by 


Tt = 


+ T. + T 
m b 


tp 


(58) 


Energy expressions in terms of generalized coordinates. - To express the strain 
and kinetic energies in terms of generalized coordinates, the displacements W and V 
are assumed to be represented by the following series: 


_ N 

V(x,t) = Yj 
1 


N M 

W(x,t) = ^ Y 

1 1 

In equations (59), represents a set of specified functions that are linearly independent 
over the range 0 x £ 1 and satisfy, at least, the geometric boundary conditions for 
the boom. Also, the time -dependent functions represent the unknown generalized 
coordinates. Note that equations (59) satisfy the compatibility of displacement at the tip. 
Usii^ the assumed solution given by equations (59) in the expression given for strain 
energy by equation (54) yields, after evaluating the resulting integrals, the following 
equation for strain energy: 



where and 6^^ (Kronecker delta function) are defined as follows: 
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(61) 


^mn = Vn dx - (l + 


0 for m = n 


mn 


[i - (-1)“™] 


V 


2 

2 

m 

+ n 

f 2 

2\2 


-n jj 


(62) 


for m # n 


^0 for m n 

6 =<; 
mn 

^1 for m = n 

Similarly, from equations (59) and (58) the expression for kinetic energy becomes 


(63) 


T = i M. Z ^ 
2 ^ 


N N 




1 1 


m 




n=l m=l 


_ M 

— 2^%.n| 
1 


(64) 


where 


Gmn = Vm * (j“m * %)’'n»)V<« 


(65) 


Equations (60) and (64) represent the strain and kinetic energies in terms of N + M 
generalized coordinates q^^. The equations of motion can now be written in terms of 
the generalized coordinates by a direct application of Lagrange's equation. 

Equations of motion. - Lagrange's equation is given by 

A/^\ + iIL=o for r = 1,2, . . ., N,N+1,N+2,N+M (66) 

dt yaq^y aq^ 

From equations (60) and (64), equation (66) yields for 1 £ r ^ N the following set of 
equations : 
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El 



‘iN+mr"® r = l,2,. . .,N (67) 


Similarly, for N + l:Sr^N + M these equations yield 



Equations (67) and (68) represent the equations of motion in terms of generalized coor- 
dinates. In matrix notation these eqiiations can be expressed as 


^Mfq} + [K]{q} = {0} 


(69) 


where [m] and [ic] are (N + M) x (N + M) symmetric, nondime nsional mass and stiff- 
ness matrices, respectively, and {q} is a (N + M) x l vector of the generalized coor- 
dinates . Equation (69) can be expressed in partitioned form as 
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1 — 12 


'‘^2 


■ 

1 k12‘ 

(N xN) 

' (N X M) 

i Sj22 
i 


5n_^ 

^N+l 

K + 

(N xN) 

“i2"r 

1 

1 (N X M) 

I |j22 

1 

(M xN) 

1 (M X M) 


• 


(M xN) 

1 (M X M) 
1 




^M+N 



1 j 




0 

^2 


0 



0 

- y 

%+l 


0 

%+N 


0 



J 


( 70 ) 


where the elements of the stibmatrices and m, n = 1, 2 are given, using 

equations (67) and (68), by the following: 



for i, j = 1,2, . . .,N 


(71a) 

IJ 2 

for i, j = 1,2, . . . ,M 


(71b) 

i i ^ j 

for i = 1,2,. . , ,N; j = 1,2,. . 

. ,M 

(71c) 

k11=k.. 
1] 1] 

for i,j = 1,2, . . .,N 


(71d) 


it2 

+ — b.. for i, j = 1,2, . . . ,M 

P 


(71e) 


1 for i = 1,2,. . .,N; j = 1,2,. 

. .,M 

(71f) 


When the set of N assumed modes of the boom i7j^(x) are specified, equation (71) can be 
used to generate the mass and stiffness matrices for any value of M. Equation (69) 
represents a set of linear, second -order differential equations whose solution is given by 

{q} = -w2{q} (72) 
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where o> is the angular natural frequency of the system. Substituting equation (72) into 
equation (69) and multiplying throi^h by (k]“^ yield 


[K]"^[M]{q} =-^{q} 


(73) 


Equation (73) represents an eigenvalue problem whose solution is discussed in numerous 
texts on vibration theory (see e.g., ref. 9 or ref. 13). 

Mass and stiffness matrices using cantilevered beam mo de sh apes. - As previously 
stated, the success of equation (73) in predicting the mode shapes and frequencies of the 
array depends primarily on the choice of assumed modes f?j^(x) . An obviously good 
choice for this generating set that satisfies all the required geometric boundary condi- 
tions would be the modal functions of a uniform cantilevered beam. In this section the 
mass and stiffness matrices given by equation (70) will be developed based on these as- 
sumed modes. The mode shapes for a uniform cantilevered beam are given by refer- 
ence 14 

^^^^(x) - cosh(^pj^x) - cos(^j.j^x) - 

where /3 is the frequency parameter for a cantilevered beam and is a function 

^cn' “^cn angular natural frequency of the beam, 



M, 9 
El 




sinh - sin <3^^ 

cosh + cos 


(75) 


Values for and as given in reference 14 are included in this report in table I. 

From equation (74) and the well-known orthogonality conditions for uniform beams, the 
following relations exist: 



Vm = °mn 







n"7?” 

'n 'm 


dx = 


^mn^cn J 


( 76 ) 
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It can also be shown after a considerable amount of manipulation that the following rela- 
tions are valid for the assumed mode shape: 


ry^(l) = 2(-l)"+l 



From equations (61), (65) ^ (71), (76), and (77) the elements of the mass and stiffness 
submatrices are given by the following: 


M?-.^ = 6.. + 4/1m + M, 

1 ] 1 ] m ' 




1] 


^ j 


K,V = (5ei - ^ (l * ^)j 


ic” = 

1] 




V Cl ^ ’ f^C] 


k}? = m [l .(.l)j]l 


1] ]i 


nP 3 


(78a) 

(78b) 

(78c) 

(i?^ j) (78d) 

(78e) 


From equations (78), (71b), and (71e) the mass and stiffness matrices can be determined 
when the mass ratios and and the load parameters k and P are specified . 

Note that these matrices are valid for both the 1-g and 0-g configurations. For the 0-g 
case, 1/P is set equal to zero in equations (78c), (78d), and (78e) . 

Two-mode approximation using preloaded cantilevered beam mode . - For the 0-g 
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configuration a good approximation for the first two lowest bending frequencies can be 
obtained by using an assvuned mode corresponding to the first mode of a uniform canti- 
levered beam under directed axial load (ref. 8). It can be expected that this type of as- 
sumed mode will give more accurate res\ilts than the results of the previous section for 
a two-mode approximation, since it more nearly approximates the condition of the boom 
in the array. Now for N = 1 and M = 1 the solution to equation (70) becomes 


(Kj, - ^T?j(l)q2 = 0 

— T7j(l)qj + i = 0 

7T 2 




> 




(79) 


Equating the determinant of the coefficients of the unknowns in equations (79) to zero 
yields the characteristic equation given by 

+ C 2 = 0 (80) 


where 


^2 = - 2 


M _ 
— ^l(l) 

7T 


(81a) 


Bg = (81b) 

C 2 = TT^Kjjk^ (81c) 

Equation (80) is a quadratic equation in that can be solved once the coefficients A 2 , 
B 2 , and C 2 are known. These coefficients, which are given by equations (81), can be 
written in terms of the beam mode by using equations (61) and (65) for Gjj and Kjj. 
From equation (65) 


Gii- ( 82 ) 

Next from equation (61) and the orthogonality relations for a cantilevered beam under 
directed axial load. 
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K 


11 


= Pdi C ”1 


(83) 


where is the frequency parameter. Substituting equations (82) and (83) into equa- 


^ 2 

tions (81) and noting that k= ^ 0 


^1 


Bg - 



M /i - + M 

ml 3 2 ) tp 

5?(i) 

\ 7T / 

J 


/ 


cr ,0 




dx + 77 '^ ^ 


C2 = 


cr ,0 




7 ]\ dx 


(84a) 


(84b) 


(84c) 


Equations (84) describe the coefficients in equation (80) in terms of the parameters 

^1 _g 

^dl’ ^ ^ cantilevered beam under directed axial load. As shown 

in reference 8 these parameters are a function of k or P/P^j. q- Table II, which was 
developed from the results of reference 8 , gives the values for these parameters as a 
function of P/P^j. Q' from this table the coefficients A 2 , B 2 , and can be 

determined for ttie various values of P/P^j. q’ first and second lowest fre- 

quencies can subsequently be determined from equation (80) . 


RESULTS AND DISCUSSION 
Continuum Mechanics Approach 

To determine the effect of the various parameters on the frequencies of a split- 
blanket solar array in both a 0-g and a 1-g cantilevered configuration, equations (33), 
(38), (42), and (43) were solved numerically by the method of bisection (ref. 15). The 
results for the first three bending frequencies for the 0-g case are shown in figures 5, 

6 , and 7. These figures show the variation in frequencies as a function of the axial load 
ratio P/Pgj. Q and the mass parameters and The tip-piece mass param- 
eter Mj.p was varied between 0 and 40 and the blanket mass parameter was varied 

between 1 and 6 . This should cover the useful range of these parameters. The axial 
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load ratio P/P„„ « represents the ratio of boom load to buckling load of the array in a 

Ci j 2 2 

0- g configuration; that is, 0 “ ^ when P/Pgj. q “ ® (blanket tension 

of zero) and when P/P^j. q = the frequency parameter jSj is zero. It follows then 
that at some point between 0 and 1, reaches a maximum. This implies that for a 
given array (M^ and fixed) there is an optimum preload P that will yield a max- 

imum fundamental frequency. This, of course, may not correspond to the minimum 
frequency of the array since the fundamental torsional frequency must also be con- 
sidered. However, it does appear that this optimum preload may be significant in de- 
signing a lightweight solar array . 

Similarly, figures 8 and 9 illustrate the variation in the fundamental and second 
bending frequencies of a solar array suspended vertically in a 1-g field. In addition to 
the parameters used for the 0-g case, an additional parameter must be specified for the 

1- g configuration, namely, P^^ O^^m’ graphs shown in these figures. 

Per o^m ~ figure 8 that when P/P^j. q ^ "^his implies 

that ’the buckling load of the array is given by P^^. « i. 52 P^^ q. This result has been 
checked by solving the buckling equation for the array given in ’appendix A by equa- 
tion (A14) . This equation has also been solved for various values of W^/P^^ q, and 
the results are presented in figure 10. This figure gives the buckling load of an array 
suspended vertically in a 1-g field as a function of blanket weight. 

Figures 11 and 12 show the effects of both tip-piece mass inertia and torsional beam 
stiffness on the fundamental and second torsional frequencies of an array in a 0-g and a 
1-g configuration, respectively. As in the bending case, the torsional frequencies of an 
array in a 1-g field are illustrated by taking P^.^. o/^m ~ ^ • 

Finally, figure 13 illustrates the effect of axial load on the mode shapes of a solar 
array in a 0-g field. The mode shapes were determined from equations (44) and (48) for 
values of = 3 and = 1 . These figures show that the mode shapes can be very 
dependent on the axial load. For small loads (figs. 13(a) to (c)), the first three modes 
are predominately blanket modes. For the point where the fundamental bending fre- 
quency is a maximum (figs. 13(d) to (f)), both blanket and boom participate equally in the 
modal displacement. For larger loads (figs. 13(g) to (i)), the first mode is predomi- 
nately a beam mode, but the second and third modes show a coupling effect. 


Rayleigh -Ritz Approach 

The Rayleigh -Ritz method presents an alternate and somewhat simpler approach for 
determinii^ the modes and frequencies of a split -blanket solar array. The accuracy of 
this method depends primarily on the type and number of assumed functions used to gen- 
erate the required mass and stiffness matrices given in equation (70) . The accimacy of 
this method, which uses cantilevered beam modes for the boom and strii^ modes for the 



24 



blanket, was investigated by solving the eigenvalue problem given in equation (73). This 
eqiiation was solved for various mass ratios and and axial load ratios 

P/Pp^ Q. Both 0-g and 1-g configurations were investigated, 

Tlie eigenvalue problem was solved by matrix iteration to find the lowest, or domi- 
nant, frequency parameter of the system. The higher modes and frequencies were de- 
termined by the technique known as "deflation" of a matrix. This method uses the 
orthogonality conditions that exist between the modes to sweep out the known, lower 
modes from the dynamic matrix given by K“^M. This technique is well described in 
references 9 and 13. 

Typical results and comparisons with exact results for the 0-g and 1-g configura- 
tions are shown in tables III and IV. These tables indicate that good results, at least for 
the first three frequencies, can be obtained in most cases by taking three beam modes 
for the boom and three string modes for the blanket. These tables also illustrate the 
improvement in accuracy achieved by increasing the number of assumed modes. 

In addition to the solutions for equation (73), table III also gives the solution to 
equation (80) . This equation is based on a two-mode approximation that uses the first 
mode of a cantilevered beam under a directed axial load and the first string mode. Pa- 
rameters required in the solution of equation (80) relating to the modal characteristics 
of a beam under directed axial loads were developed from the equations presented in 
reference 8 and are given in table III, The exact results are compared with the two- 
mode approximate results in table III and there is good agreement, for the fundamental 
frequency, for all values of q. Also, for the second frequency the error is less 

than 1 percent, provided that P/Pg^. q > 0.08. 

CONCLUDING REMARKS 

Two methods have been presented for calculating the modes and frequencies of a 
large split -blanket solar array in both a 0-g and a 1-g cantilevered configuration. The 
first method is based on the equations of continuum mechanics; the second method is 
based on the Rayleigh -Ritz, or assumed mode, approach. 

The continuum mechanics approach results in a highly transcendental equation that 
must be solved by numerical techniques. This has been done for a wide range of solar 
array parameters, and the results have been presented in graphical form. 

The Rayleigh-Ritz approach uses cantilevered beam modes to represent the boom 
displacement and string modes to represent the blanket displacement. Based on these 
assumed modes the mass and stiffness matrices of the array were developed for an 
arbitrary number of assumed modes. The resulting eigenvalue problem was then solved 
by matrix iteration for an increasing number of assumed modes. The results were then 
compared with the exact results of the first method. This comparison disclosed that 
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good results are obtained for the lower bending frequencies by using a relatively small 
number of assumed modes. 

Lastly, a two-mode approximation for an array in a 0-g configuration was pre- 
sented. The assumed modes used for this case consisted of the first cantilevered mode 
for a uniform beam under directed axial load and the first string modes. The charac- 
teristic equation for this case is a quadratic equation that can be solved in closed form 
for the first and second bending frequencies. These results compared well with the 
exact solution for all the cases investigated. This type of solution can be used to deter- 
mine the first two lowest bending frequencies of a solar array very rapidly without the 
aid of a computer . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 22, 1976, 

506-22. 
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APPENDIX A 

BUCKLING LOAD FOR A SOLAR ARRAY IN A 1-g FIELD 

The buckling load for a solar array in a 1-g field can be determined from the equa- 
tions of motion with the acceleration terms set equal to zero. From equations (2), (6b), 
(7), and (11) these equations in nondime nsional form are 



±(t^\=o 

(Al) 


dx V dx / 




(A2) 


dx^ dx^ 


d^7(l) 

, P2 d?(l) p d¥(l) _ Q 

(A3) 

dx^ 

dx dx 


where k = PZ^/EI and T is the nondimensional blanket tension given by equation (4) as 


_ T _ 

T = -^ = 1 + 1 - £ (A4) 

P p p 


The boundary conditions and compatibility relation for this case are 

w(0) = 0 

7(0) =0 -) 


dv(0) _ 0 
die 




- 0 

d^2 


w(l) = v(l) 


(A 5) 


(A6) 


(A7) 


The buckling load can be determined by solving equations (Al) and (A2) subject to condi- 
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tions (A3) and (A5) to (A7) . Thus, from equations (Al) and (A4) 


dw _ C 

dx 1 + i _ iL 
P P 


(A8) 


Integrating this equation and applying the boundary condition given by equation (A5) 


■w(x) = CP In 



(A9) 


Next the solution to equation (A2) is 


v(x) = A sin kx + B cos kx + Dx + E 


(AlO) 


Using the boundary conditions given by equations (A6) gives the following relations; 


A = -E 


B = ^ tan k >■ 
k ' 


E = - E tan k 
k ^ 


(All) 


Substituting equations (All) into (AlO) yields 


v(x) = E [tan k(cos k - 1) - sin kx + kx] (A 12) 

k 

Equations (A9) and (A12) now contain two unknowns, C and D, which can be evaluated by 
using the two remaining conditions given by equations (A3) and (A 7) . The condition 
given by equation (A3) yields 


C = D (A13) 

and using (A13) in addition to (A9) and (A12) in equation (A7) yields 
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(A14) 


tan k = k jl - P In ^1 + 


Equation (A14) is the characteristic equation from which the critical buckling load P^^ 
can be determined. Note that when P— °° (i.e., 0) the right side of equation 

(A14) becomes zero and the solution to (A14) yields 


cr,0 


jt^EI 


(A15) 


which is the critical buckling load of the array in a 0-g field. Expression k and P in 
terms of P^^ q gives the following relations : 


P = -Z. 


P \rcr_J3^ 




— Pi' 
k = — 


El 


cr,oy 


^ Pcr.O^' 

El 


= 7T 


cr,0 


(A16) 


By using (A16) in (A14), q can be determined as a function of W^/P^^. q. 
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APPENDIX B 


ASYMPTOTIC FORMULAS FOR Fj(/3), F2(^), Fg(^p, AND 

The characteristic equations for bending and torsion of a solar array in a 1-g field 
contain certain functions that can be evaluated more simply for large values of their 
arguments. This will be the case when the blanket weight approaches zero. To 

simplify the writing of these equations, let 


Zj = 


7t2 


where 


(Bl) 


M^/^P 
^ _ m 


(B2) 


In terms of Zj, equations (30) and (32) become 

Fj(J) = Jq(2Zj)Yq(2Zj?q) - Yq(2zPJq(2Zj?q) 


^ 7T J 


(B3) 


(B4) 


Now, as approaches 0, approaches °° and the arguments of the Bessel func- 

tions in equations (B3) and (B4) become large. Thus, the asymptotic formulas for these 
functions can be used. For large values of Zj these formulas as given by reference 10 
are 


,(2z^)=^cos(2z,-l) 


Yq(2Zj) = Jj(2zj) 


(B5) 


Similar expressions exist for Jq(2z^Cq) and Y(2 Zj?q). From equations (B 5) and (B 3), 
Fj(jS) reduces to 
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Now for small values of or large values of P, equation (B2) yields for the term in 

brackets in equation (B6) 


~ (^0 ” 


2 - 




w a. 


where 

m 


From equations (B6) and (B7) 


(B7) 


(B8) 






Sin a , 


(B9) 


Next consider equation (B4) for large values of Zj. From equations (B5) this equation 
becomes, after a little manipulation, 


FgO) = ^cos[2Zj(?q - 1)] (BIO) 

—2 

The terms in brackets is again given by equation (B7) . Noting that /Cj^/3 = cxgP (from 
eqs. (B2) and (B8)), equation (BIO) can be written as 




Q-gP 




cos a. 


(BID 


In a similar straightforward manner it can also be shown that Fg(/3j.) and F^(jSj.) as given 
by equations (36) become for large values of P 



I Hill mil II I I 


Fg(j3t) - 


37/C 

37/C2^0 


• Sin 





(B12) 


32 



APPENDIX C 


ORTHOGONALITY CONDITIONS 

The orthogonality relations for the bending modes of a solar array in a 1 -g field can 
be developed by starting with equations (18) and (19) and the corresponding boimdary con- 
ditions given by eqxiations (20) and (21) . For an array vibrating in its i*^*^ mode these 
equations become 



(Cl) 


dVj _ 4 _ 

L + 1 


dx^ 


dx^ 


= 0 


(C2) 


d^v.(l) 

d5^ 




dv^(l) 

dx 


2 P d? 




= 0 


(C3) 


with boundary and compatibility conditions 


v.(0) - 0 

dv.(O) 
_^= 0 
dx 


(C4a) 

(C4b) 


Wi(?o) = 0 
w.(l) = v.(l) 


dV^(l) 

d5i^ 


(C4c) 

(C4d) 


(C4e) 


First operate on equation (Cl) by multiplying it by Wj, where j corresponds to the j 
mode of vibration, and integrate over the length. This yields, after an integration by 
parts and applying conditions given by (C4c) and (C4d) , 
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(C5) 


Vj(l) 


dw.(l) 

dC 



dw. dw. _ 2 4 

— i. — _1 dx + 

d? d? ^ ^ 



?WjWj dx = 0 


Next, multiply equation (C2) by Vj and integrate over the length to obtain 



Integrating the first two integrals by parts and applying the boundary conditions given by 
equations (C4a) to (C4c) yield 




The final step is to multiply equation (C5) by ^/2P and then add the resulting equation 
to equation (C7) . The result, after noting that k^/Cj/^ = Mj^P, is 
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If i and j are interchained in equation (C8) , a second equation is obtained that is simi- 
lar to (C8) except that ^ is replaced by Thus, 



Subtracting equation (C9) from (C8) and defining the generalized mass as 


M. = ^ dx + 2M^P d? + M^pvf (1) (CIO) 

yields 

/'v,Vj<h.2M„P A ?WjW. dx + Mj.pVj(l)v.(l) = 6„M. (Cll) 

Equation (Cll) is the first orthogonality condition for the bending modes of a solar array 
in a 1-g field. A second orthogonality condition can be obtained by dividing equations 
(CIO) and (C9) by ^ and respectively, and subtracting the resulting equations. 

The result is 
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dVj dVi 

(S’ dx 


Iv •‘■J ■*• 1 


(C12) 


To obtain the orthogonality conditions for the array in a 0-g field, ? is transformed to 
X by using equation (4) and the limit is taken as approaches zero. The results 
yield the first and second orthogonality conditions in the following form: 


VjVj dx . WjWj dx + MjpVj(l)Vjd) = (C13) 


/•I 

/ dw. dw. _ 

/ 

/ dx dx 


d\. d^v. _ _2 

—A dx - 
dx dx 


dx dx 


dx = 


(C14) 


where for this case is given by 


M 


f ^ dx + M f dx + M._vf (1) 
1^1 m ^ 1 tp i' ' 


(Cl 5) 


The orthogonality conditions for the torsional modes can be developed in a manner simi- 

'f'h 

lar to that used for the bending modes. For the i*^*^ mode of vibration, equations (22) 
to (24) yield the following set of equations and boundary conditions: 



J.4i!i\ + 4/32pC<p =0 

d?\ d?/ “ ^ 

(C16) 


dx"^ 

(C17) 

de^d) 

dx 

1 — -2—2 
- ±^— 2 — - L„kM0 (1) = 0 
24 P dC tp t ti 1 

(CIS) 
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<P(Cq) = 0 

0 .( 1 ) = <^.( 1 ) 


(C19) 


0(0) = 0 (C20) 

Multiplying equation (C16) by <Py integrating over the length, and applying the boundary 
conditions given by equations (C19) yield 




%/¥> 


(C21) 


Similarly, multiplying equation (C17) by 0^, integrating over the length, and using equa- 
tions (CIS) and (C20) give 


1 d«p(l) d0. d0. _2 _2 

-!-^0.(l — i / ^-^cE+(3lkf / 0.0. d5T + I. 0,(1) 0,(1) 

24 P J d? / dx dx t / 1 ] tp ^ ' 

«'0 


= 0 


(C22) 


Now multiply equation (C21) by k^/24P and subtract from equation (C22) to obtain 



Interchanging i and j in equation (C23) yields 
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Subtracting equation (24) from equation (23) yields the first orthogonality condition for 
the torsional modes of the solar array in a 1-g field. Thus, 


1 ^ 

fb / d£ * £ / 0 d? 

^0 o yj 


= 6..T. 
1] 1 


(C25) 


where is the generalized mass defined by 


T. 

1 




0^ dx + 


P 

6 / 






(C26) 


Similarly, the second orthogonality condition for this case is obtained by dividing eqxia- 
tions (C24) and (C23) by j3y and respectively, and subtracting the resulting eqxia- 
tions. The results can be expressed as 



d0. d0. 4 

_^^dx + -4 

dx dx 24P 



dcp. d(p. _ _ 2_2 

?^_ldx = 6 kf/3fT 
d? d^ 1] t ti 1 


(C27) 


Finally, the orthogonality relations for the 0-g configuration can be obtained by equa- 
tions (C25) and (C26) by using equation (4) to transform ^ into x and taking the limit 
as approaches zero (i.e. , P —°o) . The results are 

Tb / ejO. dx + ± dx . T,pe,(i) 0.(1) = 6..T. (C28) 
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(C29) 


where 



d9. d9. 
dx dx 



d(p. dcp. _ _2_2 

— I — 1 dx = 6,,k;i3f.T. 
dx dx ij t ti 1 


in equations (C28) and (C29) is given by 


T. 

1 



0? dx + 


12 
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APPENDIX D 


SYMBOLS 


^ 2 ’ ® 2’^2 
b 

b 

b 

mn 

El 


G 


mn 


y 

% 

JG 

[K] 

k 

K 


^b 

M. 

1 


M 


m 


M. 

tp 

M^(x) 

My(x) 


M. 


m 


tp 


arbitrary constants of integration 
constants defined by eq. (81) or (84) 
blanket width 

nondimensional blanket width, h/l 

elements of stiffness matrix defined by eq. (62) 

boom bending stiffness 

element of mass matrix defined by eq. (65) 

mass polar moment of inertia per unit length of boom 

mass moment of inertia of tip piece about its center of gravity 

2 

boom inertia ratio, 

2 

tip -piece inertia ratio, 
boom torsional stiffness 

element of stiffness matrix defined by eq. (61) 

nondimensional stiffness matrix 

i^, 3^^ submatrix of the stiffness matrix K 

axial load parameter for bending, Vp^^/ei 

torsional stiffness factor, "^jpb^/JG 

blanket and boom length 

total mass of boom, p^l 

i*^*^ generalized mass of solar array for bending 
total blanket mass, p^l 
mass of tip piece 

torsional moment distribution along boom 
bending moment distribution along boom 
mass ratio, 
mass ratio, 
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[m] 

[Spj] 

P 

P 

Per 

Pcr,0 

Q(x) 

q„(t) 

{q] 



m 

V(x,y,t),v(x) 

V,v(x) 


W(x,y,t),w(x) 

W,w(x) 

X 

X 


y 

cn 


nondime nsional mass matrix 

i submatrix of mass matrix M 

uniform stress resultant in blanket at rest 

compressive preload in boom 

ratio of axial load to blanket weight, P/Wj^ 

critical buckling load of array 

critical buckling load of array in 0-g configuration 

shear distribution along boom 
th 

n ^ generalized coordinate 

column matrix of generalized coordinates 

blanket tension per unit width 

kinetic energy function of boom 
th 

i generalized mass of solar array for torsion 

kinetic energy function of blanket 

total kinetic energy of array 

kinetic energy of tip piece 

total strain of array 

strain energy function for boom 

strain energy function for blanket 

boom displacement 

nondimensional boom displacement, V/Z and v/l 
total blanket weight, 
blanket displacement 

nondimensional blanket displacement, W/Z and w/Z 

longitudinal coordinate of blanket and boom 

nondimensional longitudinal coordinate 

lateral coordinate of blanket 
th 

n modal parameter for uniform cantilevered beam 
characteristic values, eqs. (27) and (41) 
bending frequency parameter, ^(M^,Z VeI)o;2 



IIHH I L JH II I I IMII 


/3gj^ natural frequency for uniform cantilevered beam 

first natural bending frequency parameter for uniform beam under directed 
axial load 
til 

|3jj n natural bending frequency parameter 

/3^ torsional frequency parameter 

/3^ n*"^ natural torsional frequency parameter 

Kronecker delta function 

C transformed coordinate defined by eq. (4) 

?Q transformed coordinate at point x = 0 

i?j^(x) n'"” assumed function for boom 

0(x,t), 9{x) rotation of boom cross section 

bending parameter, P/k^M^ 

«2 torsion parameter, jS^'P 

mass per unit length of boom 

mass per unit length of blanket 

<P{^) torsional displacement function 

u) circular frequency of vibration 

th 

n*^*^ natural circular frequency of uniform cantilevered beam 

first natural frequency of uniform cantilevered beam under directed axial 
load 
th 

Wjj n'^ natural circular frequency of vibration 

Subscripts: 


differentiation with respect to x 
differentiation with respect to time t 
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TABLE n. -STRUCTURAL DYNAMIC CHARACTERISTICS 


TABLE I. - VALUES OF W AND 

^cn cn 


FROM REFERENCE 14 


n 

th 

n natTiral 
torsional fre- 
quency para- 
meter, 

^cn 

modal pa- 
rameter for 
uniform beam, 

^cn 

Remarks 

1 

1.87510407 

C 

). 7340955 

For n > 5 

2 

4.69509113 


1.01846644 


3 

7.85475743 


.99922450 

«n 1 

! 4 

10.99554074 

1.000033553 


5 

14.13716839 


.9999985501 

1 

6 


1 

1 

7 

i »(t) 



1 

i 

8 

1 



1 

9 

1 ‘’(1) 



( 

! 

10 






OF CANTILEVERED BEAM UNDER AXIAL 


LOAD DIRECTED THROUGH ROOT 


Axial load 
ratio , 

P/Pcr,0 

Tip displace - 
ment. 

Generalized 

mass, 

Frequency 

parameter, 

^dl 

0 

2.00 

1.00 

1.875104069 

.001 

2.000687704 

1 . 000542243 

1.874861498 

.01 

2.007105600 

1.005646573 

1.872665819 

.02 

2.014727517 

1.011801344 

1.870199290' 

.04 

2.031574652 

1.025702328 

1.8651787901 

.06 

2.050653505 

1.041844802 

1.860037145 

.08 

2.072081985 

1.060384811 

1.854768565 

.10 

2.095984442 

1.081494411 

1 . 849366860 

.20 

2.257519322 

1.233025493 

,1.820112101 

.30 

2.505117303 

1.488716586 

'1.7862955231 

.40 

2.870244314 

1.910797481 

,1.746376307 

.50 

3.407887443 

2.624651093 

11.697963692 

.60 

4.227974273 

3.919604287 

' 1.637089120 

.70 

5.588861642 

6.609907291 

1.556494718 

.80 

8.269873146 

13.874364178 

1.440685832 

.90 

16.176957259 

50.465329261 

1.244925094 

.93 

22.906582764 

99.461655188 

11.148904957 j 

.96 

39.687051208 

293.170363886 

1.008218686 

.98 

78.7905575842 

1140.85733402 

.853258066 

1.00 

OO 

OO 

0 1 
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TABLE m. - CX)MP ARISON BETWEEN EXACT SOLUTION AND RAYLEIGH -RITZ 


SOLUTION FOR BENDING FREQUENCIES OF SOLAR ARRAY IN 0-g FIELD 

{Ratio of blanket mass to total boom mass, M^, 3.0; ratio of tip-piece mass to 
total boom mass, 1.0.] 

(a) Fundamental bending frequency 


\xial load 
ratio, 

^/^cr,0 

Exact- 

solution 

equation 

Rayleigh-Ritz solution using cantilevered 
beam modes 

Two -mode 
Rayleigh-Ritz 
solution 
(eq. (80)) 


Error, 

per- 

cent 

N=4, 

M=4 

! 1 

Error, 

per- 

cent 

N=5, 

M=5 

Error, 

per- 

cent 

! 

N=l, 
j M=1 

Error, i 
per- 
cent 

0.001 

0.42378 

0.42378 

0 

' 0.42378 

0 

0.42378 

0 

0.42380 

0.005 

.01 

.74051 

.74054 

.004 

.74053 

.003 

.74052 

.001 

.74105 

.073 

.02 

,85891 

.85906 

.012 

1 .85901 

.006 

.85899 

.003 

.86063 

1 .195 

.04 

.96058 

.96085 

.028 

.96070 

.012 

1 .96064 

.006 

.96495 

.455 i 

.06 

, .99854 

.99893 

; .039 

.99871 

.017 

i .99863 

.009 

1.0045 

.597 

.08 

1.0149 

! 1.0153 

.039 

1.0151 

.020 

! 1.0150 

.010 

1.0214 

.640 

.10 

1.0224 

' 1.0229 

.049 

1.0226 

.020 

' 1.0225 

■ .010 

1.0291 

1 .655 

.20 

,1.0225 

: 1.0231 

.059 

1.0227 

.4)20 

j 1.0226 

.010 

1.0289 

.626 

.30 

' 1.0052 

I 1.00607 

1 .087 

1.0055 

.030 

1 1.0054 

.020 

' 1.0112 

.597 

.40 

.98085 

.98206 

.123 

.98125 

.041 

.98110 

.025 

.98631 

.557 

.50 

.95009 

i .95179 

1 .179 

' .95066 

1 .060 

.95044 

t .037 

1 .95496 

' .513 1 

.60 

.91133 

; .91380 

.271 

.91216 

! .091 

.91184 

.056 

.91550 

.458 

.70 

.86086 

.86462 

.437 

.86217 

.152 

,86166 

.093 

.86419 

.387 

.80 

.79049 

.79670 

.786 

.79274 

.285 

' .79183 

.170 

.79280 

.292 

.90 

.67648 

.68909 

1.86 

.68132 

.715 

.67929 

.415 

.67760 

.166 

.93 

.62223 

.63954 

2.78 

.62905 

1.10 

.62616 

.632 

.62298 

.121 

.96 

.54412 

.57128 

4.99 

. 55520 

2.04 

. 55049 

1.17 

. 54451 

.072 

.98 

.45935 

. 50382 

9.68 

.47842 

4.15 

. 47038 

2.40 

.45952 

.037 



(b) Second bending frequency 


0.001 

0.59928 

0.59928 

0 

0.59928 

0 

0.59928 

0 

1.1730 

95.7 

.01 

1.0349 

1.0353 

.039 

1.0351 

.019 

1.0350 

.010 

1.1912 

15.1 

.02 

1.1559 

1.1573 

.121 

1.1565 

.052 

1.1562 

.026 

1.2178 

5.35 

.04 

1.2606 

1.2619 

.103 

1.2612 

.048 

1.2609 

.024 

1.2873 

2.12 

.06 

1.3463 

1 . 3472 

.067 

1.3467 

.030 

1.3465 

.015 

1.3639 

1.31 

.08 

1.4225 

1.4233 

.056 

1.4228 

.021 

1.4227 

.014 

1.4363 

.970 

.10 

1.4900 

1.4907 

.047 

1.4903 

.020 

1.4902 

.013 

1.5019 

.799 

.20 

1.7424 

1.7430 

.034 

1.7427 

.017 

1.7426 

.011 

1.7516 

.528 

.30 

1.9188 

1.9194 

.031 

1.9191 

.016 

1.9190 

.010 

1.9278 

.469 

.40 

2.0571 

2.0576 

.024 

2.0573 

.010 

2.0572 

.005 

2.0664 

.452 

.50 

2.1721 

2.1727 

i .028 

2.1724 

.014 

2.1722 

.005 

2.1819 

.451 

.60 

2.2713 

2.2719 

.026 

2.2716 

.013 

2.2715 

.009 

2.2817 

.458 

.70 

2.3590 

2.3596 

.025 

2.3593 

.013 

2.3591 

.004 

2.3702 

.475 

.80 

2.4378 

2.4385 

.029 

2.4381 

.012 

2.4380 

.008 

2.4500 

.500 

.90 

2.5097 

2.5103 

.024 

2.5100 

.012 

2.5098 

.004 

2.5229 

.526 

.93 

2.5300 

2.5307 

.028 

2.5304 

.016 

2.5302 

.008 

2.5436 

.538 

.96 

2.5499 

2.5506 

.027 

2,5503 

.016 

2.5501 

.008 

2.5639 

.549 

.98 

2,5629 

2.5636 

.027 

2.5633 

.016 

2.5631 

.008 

2.5772 

.558 


(c) Third bending frequency 


0.001 

0.73387 

0.73388 

0.001 

0.73387 

0 

0.73387 

0 



.01 

1.2127 

1.2177 

.412 

1.2146 

.157 

1.2137 

.082 



.02 

1.3489 

1.3517 

.208 

1.3501 

.089 

1.3495 

.044 



.04 

1.5681 

1.5696 

.096 

1.5687 

.038 

1.5684 

.019 



.06 

1.7267 

1.7279 

.069 

1.7272 

,029 

1.7270 

.017 



.08 

1.8515 

1.8526 

.059 

1.8520 

.027 

1.8518 

.016 



.10 

1.9554 

1.9565 

.056 

1.9559 

.026 

1.9556 

.010 



.20 

2.3202 

2.3213 

.047 

2.3206 

.017 

2.3204 

.009 



.30 

2.5656 

2.5668 

,047 

2.5661 

.019 

2.5659 

.012 



.40 

2.7555 

2.7568 

.047 

2.7560 

.018 

2.7558 

.011 



.50 

2.9122 

2.9136 

.048 

2.9128 

.021 

2.9125 




.60 

3.0462 

3.0477 

.049 

3.0468 

.020 ! 

3.0465 




.70 

3.1630 

3.1647 

.054 

3.1638 

.025 

3.1634 

.013 



.80 

3.2643 

3.2665 

.067 

3.2654 

.034 

3.2648 

.015 



.90 

3.3411 

3.3464 

.159 

3.3441 

.090 

3.3424 

.039 



.93 

3.3525 

3.3606 

.242 

3.3573 

.143 

3.3547 

.066 



.96 

3.3542 

3.3661 

.355 

3.3613 

.212 S 

3.3573 

.092 



.98 

3.3498 

3.3641 

.427 

3.3583 

.254 

3.3535 
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TABLE IV. - COMPARISON BETWEEN EXACT SOLUTION AND 


RAYLEIGH -RITZ SOLUTION FOR BENDING FREQUENaES 

OF SOLAR ARRAY IN A 1-g FIELD 

[Ratio of blanket mass to total boom mass, M^, 3.0; ratio of tip- 
piece mass to total boom mass, 1-0; ratio of array 0-g 

critical buckling load to total blanket weight, O^m’ 

(a) Fimdamental bending frequency 


Axial -load 
ratio, 

^^cr,0 

Exact- 

solution 

eqviation 

I Rayleigh -Ritz solution 

N=3, 

M=3 

Error, 

per- 

cent 

N=5, 

M=5 

Error, 

per- 

cent 

N=10, 

M=10 

Error, 

per- 

cent 

0.001 

1.3261 

1.3445 

1.387 

1.3416 

1.169 

1.3379 

0.890 

.050 

1.3433 

1.3453 

.149 

1.3442 

.067 

1.3436 

.022 

.10 

1.3431 

1 . 3439 

.060 

1.3433 

.015 

1.3431 

0 

.20 

1.3376 

1.3380 

.030 

1.3378 

.015 

1.3377 

.007 

.40 

1.3196 

1.3201 

.038 

1.3197 

.008 

1.3196 

0 

.60 

1.2944 

1.2952 

.062 

1.2946 

.015 

1.2944 

0 

.80 

1 . 2599 

1.2613 

.111 

1.2602 

.024 

1.2599 

0 

1.00 

1 . 2096 

1.2124 

.231 

1.2103 

.058 

1.2097 

.008 

1.20 

1.1264 

1.1327 

.559 

j 1.1279 

.133 

1,1266 

.018 

1.40 

.94022 

.96197 

2.313 

.94584 

.598 

.94102 

.085 

1.517 

.35994 

.61340 

70.4 

.47149 

31.0 

.38301 

6.41 


(b) Second bending frequency 


0.001 

1.6858 

1,9045 

12.97 

1.8534 

9.94 

1.8002 

6.79 

.050 

1.9662 

2.0073 

2.09 

1.9840 

.905 

1.9706 

.224 

.10 

2.0755 

2.0926 

.824 

2.0810 

.265 

2.0764 

.043 

.20 

2.2236 

2.2281 

.202 

2.2246 

.045 

2.2237 

.004 

.40 

2.4263 

2.4269 

.025 

2.4264 

.004 

2.4263 

0 

.60 

2.5776 

2.5778 

.008 

2.5776 

0 

2.5776 

0 

.80 

2.7021 

2.7023 

.007 

2 . 7022 

.004 

2.7021 

0 

1.00 

2.8088 

2.8093 

.018 

2.8090 

.007 

2.8089 

.004 

1.20 

2.9008 

2.9020 

.041 

2.9012 

.014 

2.9009 

.003 

1.40 

2.9349 

2.9524 

.596 

2.9400 

.174 

2.9357 

.027 

1.517 

2.8131 

; 2 . 8477 

1 1.230 

2.8225 

1 .334 

2.8146 

.053 


(c) Third bending frequency 


0.001 

2.3892 

2.7145 

13.6 

2.6037 

8.98 

2.5172 

5.36 

.050 

2 . 6878 

2.8050 

4.36 

2.7305 

1.59 

2.6964 

.320 

.10 

2.8219 

2.8896 

2.40 

2.8390 

.606 

2.8240 

.074 

.20 

3.0104 

3.0390 

.950 

3.0145 

.136 

3.0107 

.010 

.40 

3.2729 

3.2803 

.226 

3.2733 

.012 

3.2729 

0 

.60 

3 . 4606 

3.4649 

.124 

3.4612 

.017 

3.4607 

,003 

.80 

3.5027 

3.5190 

.465 

3.5066 

.111 

3.5033 

.017 

1.00 

3.3704 

3.3910 

.611 

3.3756 

.154 

3.3712 

.024 

1.20 

3.1983 

3.2221 

.744 

3.2045 

.194 

3,1992 

.028 

1.40 

3.0358 

3.0492 

.441 

3.0390 

.105 

3.0362 

.013 

1.517 

3.0540 

3.0556 

.052 

3.0545 

j .016 

3.0541 

.003 
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Figure 1. - Geometry and coordinate system for a large solar array with a split blanket. 



(a) Blanket element showing forces and displacements. 



(b) Beam element showing forces, moments, and displacements. 

Figure 2. - Forces acting on blanket and beam. 




the uniform load (P/b) dW/dx(I,t) is not shown, ) 



Figure 3. - Forces and displacements of tip piece. 


X 



Figure 4. - Bending and torsional modes of vibration. 



Second bending frequency parameter. 










Fundamental bending frequency parameter, 











Figure 12. - Torsional frequency as function of tip-piece Inertia for 
split-blanket solar array in 1-g field. Ratio of array 0-g critical 
buckling load to total blanket weight. Pj.j. 10. 
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(g) Fundamental bending mode; 0 * ^ 




Nondimensional longitudinal coordinate, x 
(i) Third bending mode; ?l?^^ 0 " 

Figure 13. - Concluded. 
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